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Abstract 

Let G = (V, E) be a random electronic network with the boundary ver- 
tices which is obtained by assigning a resistance of each edge in a random 
graph in G(n,p) and the voltages on the boundary vertices. In this pa- 
per, we prove that the potential distribution of all vertices of G except for 
the boundary vertices are very close to a constant with high probability 
f or p = cIsji anc i c > l. 

1 n. 



1. Introductions 

The connection between random walks and electrical networks can be 
recognized by Kakutani [ll|] in 1945. Doyle and Snell |7| in 1984 pub- 
lished an excellent book which explained the relations between random 



walks and electrical networks. Tetali [18] presented an interpretation of 



effective resistance in electrical networks in terms of random walks on 



underlying graphs. Recently, Palacios [17[ studied the hitting times on 
random walks on trees through electroic networks. For more informa- 
tion between random walks and electrical networks, the readers may be 
referred to @,Q[l2|. 



Let G = (V, E) be a connected undirected graph without loops and 
multiple edges. To make it a electrical network, we assign to each edge 
e = (i, j) e E resistance and the conductance of e = (i, j) is c^- = — . 

ij 

We define a random walk on the electrical network iV = (G, c) be a 
Markov chain with transition matrix P = (pij) given by 



Pij 
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with q = j)£E c iy If r ij * s ass ig ne d to be a unit resistance, it is just 
the simple random walk on G. 

It is well known that the potential and current distributions on elec- 
trical networks follow both Kirchhoff's and Ohm's laws. Further, it is 
proved (for example, see j3|) that the potential distribution follows a har- 
monic function with certain boundary conditions, i.e., the potential of 
each vertex except for the boundary vertices is just the weight-average of 
the potential of its neighbor's. Here boundary vertices are those vertices 
at which there are current flowing into or out of the network. Moreover, 



Lyons and Peres in 16] investigated the potential distributions on regular 
lattice. Curtis and Morrow 0] determined the distribution of resistors in 
a rectangular network by using boundary measurements. These works 
established the intimate connection between the random walks on graphs 
and electrical current networks. It is natural to ask what we can say if 
the random walks and electrical networks are considered on a random 
graph? 

The space G(n,p), is defined for < p < 1. To get a random element 
of this space, we select the edges independently, with probability p. For 
more information and background, the readers are referred to Bollobas 
book Recently, many other random graph models, such as small- 



world model 19], BA model [l[ etc., have been proposed to simulate and 
study the mass real world networks . 

Grimmett and Kesten 

?? 

considered the effective resistances for random electrical networks on 
random model G(n,p). In this paper we mainly investigate the potential 
distributions of the electrical network on random graph G G G(n,p) 
by assigning a resistance Cjj for each edge e = (i, j). By making use 
of the probabilistic interpretation of potential distribution on electrical 
networks as well as the fast mixing property (5] for random walks on a 
random graph G G G(n,p). A sequence of events £ n is said to occur with 
high probability (whp) if lim^oo Pr(£ n ) = 1. The main result of this 
paper can be stated as follows: 

Theorem 1.1. Let N = {G,c) be a random electrical network from a 
random graph G G G(n,p) model withp = 01 , a > 1 and C\ < Cij < C2 
for all e = G E, where < C\ < C2 are two positive constants. If 

T = {xi,x 2 , ■ ■ ■ ,xk} is the boundary set with boundary potential V Xk = 
Px k i < p Xk ^ 1 for k = 1, 2, • • • , K , then whp the potential distribution 



2 



satisfies 




+ 0((}nn)- 1 ) 



(1) 



for eachie V(G) \ Y. 

The rest of this paper is arranged as follow. In sections 2, we present 
some properties for both random graphs and electrical networks. In sec- 
tion 3, we will give a rigorous proof for the main theorem (Theorem ll.il) 
and apply this result to a generalized consensus model with finite leaders. 
In section 4, we will do some further discussions on the potential distri- 
butions in more general cases where G may be circles and small-world 
networks and Cjj may be i.i.d random variables for each G E. We 
note here that we say an event holds with high probability(denoted as 
whp) if it holds with probability 1 — o(l) as n — > oo. 

2. Preliminaries 

In this section, we introduce some properties of a random graph in 
G(n,p) and give the probabilistic interpretation of potential distributions 
on electrical networks. Let d(i) denote the degree of vertex i G V and 
let S(G) denote the minimum degree of G. 

Definition 1. A simple graph G = {V,E) is said to be proper if it has 
the following structural properties PI — P5. 
PI: G is connected. 

P2: Call a cycle short if its length is at most 1 °" n . The minimum 
distance between two short cycle is at least , lr l n . 

a In In n 

P3: G has at least one triangle, at least one 5-cycle and at least one 
7-cycle. 

P4: LetCi, C*2 (Ci < C2), K be positive constants. Let N = {G,c) be the 
electrical network on G with C\ < Cij < C 2 for each edge e = G E. 
For L C V, \L\ < K, denote by G' = G[V\L] = (V',EJ the subgraph of 
G induced by V \L. For S C V\L, denote by Ec(S, S) the set of edges 
of G' with one end in S and the other in S = V \ (L U S). If \S\ < |, 



P5: There exists a positive constant 5 > such that it follows for any 



then 



T l (i,j)eE Gf (S,S) °ii Ci 



(2) 



where ^ = £ (iii)eJ5 , dy 
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vertex i G V , 

5C\nn < d(i) < 4Clnn. 

Remark From the definition, it looks very rare for a graph to be 
proper. But there are much many graphs to be proper. In fact, we have 
the following result. 

Lemma 2.1. Let G be a random graph in G(n,p) with p = ^sJk and a 
constant a > 1. Then whp G G G(n,p) is proper. 

Proof. It follows from Lemma 1 in jH[ that whp PI — P3 hold. More- 
over, by Lemma 6.5.2 in jsj, whp P5 also holds. For P4, it is a simple 
generalization of Lemma 1 in Q. In fact, by Lemma 1 in [5j], we have 

\E G ,(S,S)\ 1 
d&(S) ~ 6' 

where do'{S) = J2 i€S da>(i) and dcii) is the degree of vertex i in G' . So 
we have whp, 

£(i,j)eg G >(s,g) c »i Ci\E G >{S,S)\ _Cj_ 

This completes the proof. □ 

Let iV = (G, c) be an electrical network. For any i 6 V let W^n 
denote a random walk on iV which starts at vertex i and let Wi^{t) 
denote the walk generated by the first t steps. Let X^it) be the vertex 
reached at step t and P}%(j) = Pr(Xi^(t) — j). If the random walk on 
N is irreducible and aperiodic, let tc n be the stationary distribution of 
the random walk on N. We also need the following Lemma which is a 
slight generalization of Lemma 2 in j5[ . 

Lemma 2.2. Let G = (V, E) be proper and N = (G, c) be the electrical 
network with C\ < c^j < C% for each e = G E. For a subset 

T = {xx,X2,-" ,%k} of V , let N' = (G',c) be the induced subnetwork 
obtained from G' = G[V \ T) = (V',E'). Then there exists a sufficiently 
large constant K > such that for all i,j G V and t > t Q = K lnn, 

\P?Uj)-*N>(j)\=0(n- 10 ), (3) 

i.e., the random walk on N' mixes in time 0(ln?r.). 



4 



Moreover, set <§ = {X it N(t) ^ T,0 < t < 2t } be the event that the 
random walk on N started from i do not reach the vertices in T for the 
first 2t steps. Then we have 



Pr(X i>N {2t ) = j\S) = (1 + O((lnn)- 1 ))Pr(X 4i7V ,(2t ) = j). (5) 



Proof. The Lemma and Proof are almost the same as Lemmas 2-4 in [5] , 
in which they only consider the simple random walk on G. In order to 
keep the consistency of our paper, we will simply give a sketch proof of 
Lemma 12.21 with emphasis on different places. 

By P3 and P4, the random walk on N' is irreducible and aperiodic 
and therefore it has a unit stationary distribution ttn,. By using P4 
instead of and the isoperimetric Inequality of Lemmas2-4 in [5] , it is easy 
to see that (EJ) holds. 

For 1 < k < K, let Ng'(xi~) be the neighborhood of x^ in G' and 
Nq'(K) = Ui<i<KNG'(xk). Let 5k be the minimum degree of vertices 
in Ng'(K). Fixing i 6 V, for j e V, let W^^VV^) denote the set 
of walks on N' which starts at i, ends at j, are of length 2to and which 
leave a vertex in the neighborhood NG'(xk)(NG'(K)) exactly m times. 



Let yift = \J. Vl%, = U, W^j and W = (w , w ir ■ ■ , w 2to ) E W, 



Pr{g) = l-Odlnn)- 1 ) 



(4) 



and 



Set 



Pw 



Pr(X itN (s) = w s ,s = 0,1,- ■■ ,2t ) 
Pr(X ijN ,(s) = w s , s = 0, 1, • • • , 2t ) 



Then 



1 > Pw > (1 



C 2 K 



This is because 



Pr(X i>N (s) = w s \X i:N (s - 1) = w 8 -i) 
Pr{X i>N ,{s) = w s \X i>N r(s - 1) = w s _i) 



> 




if w s ^^N G ,(K) 
if Ws _!GiV G ,(K). 



So 



Pr(<f) = Yl Pr(W hN (2t Q ) = W) 



m>0 W£W% 



= PwPr(W % M2h) = W) 



m>0 W&W* 
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>£p£(i-^r> (6) 

rn>0 



where 



P*=Y1 MWiM 2 to) = W) = Pr(W itN ,{2t )eWZ). 

Now fix j and write 

Pr{X ijN = j\£) = Pr{W hN {2t ) = W)Pr(<g)- 1 

m>0 WGW K 

= S P^r(^(2t ) = ^)Pr(^)- 1 . 



m>0 vi/PW K 



If we set 



K = Pr(W i>N ,(2to) 6 W%j) 

Pm,j 



Pr(X itN ,{2t ) = j) 
= Pr(Wi^i(2t ) leaves a vertex of Nqi(K) m times \X it ^f/(2t ) = j) 

then 



m>0 



We can get by the same method as Cooper and Frieze showed in Lemma 4 
in H that 



where 

Pr(W^,(2t ) e W^] 



Pmj = n„/^ ,ol n 0<m<2,l<jfe<tf. 

So we have 



Pr(X^,(2t ) = j) 

2K 



and 



m=0 
2A" 



E^ 1 -^ 111 ^ 1 )- ( 9 ) 



m=0 
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Now using equations (4), (5), (6), (7) and the fact 5k > SClnn from P5, 
we have 

2K 

Pr{$) > £p£(l"£xO m > (l-^) 2Jf -0((lnn)- 1 ) = l-O(Onn)- 1 ). 

Similarly we have 

Pr(X hN (2t ) = = (1 + O((lnn)- 1 ))Pr(X^ / (2t ) = j). 

□ 

Lemma 2.3. Let N = (G, c) be a connected electrical network. Let 
T = {x k , 1 < k < K} C V be the distinct boundary vertices and define 
for each i G V 

r Xk (i) = min{t\X ijN (t) = x k } l<k<K 

K 

Pi = y^;Px k Pr(r Xk (i) = min(T Xs (i), 1 < s < K)) 

k=l 

so that P Xk = p Xk for 1 < k < K. Then {Pi,i G V} is the same as the 
distribution of potentials when x k is set at p k for 1 < k < K. Especially 
when K = 2, i.e., we choose only two vertices Xi,x 2 as boundary points 
and set p Xl = l,p X2 = 0, then Pi = Pr(r Xl (i) < r X2 (i)) is the same as the 
distribution of potentials when X\ is set at 1 and x 2 at 0. 



Proof. For a electrical network with boundary set T, there are no current 
flow into or out of the network at vertices in V \ F. Assume that {Vi, i G 
V} be the potential distribution of vertices in N. By using Kirchhoff's 
current law and Ohm's law, we have for i G V \ T 



Vi-Vj 

which implies 



i.e., the potential distribution follows a harmonic function. For fixed k 
and each i G V, set 

pf = P r (r Xk (i) = mini< s < K {T Xs (i)})- 



7 



Then 



P k = 1 



P x * = for x s e r, s ^ k. 

Moreover, if we consider the very first step of the random walk on N 
started at i G V \ V, then 



p i = E = E f p r 

(i,j)€E {i,j)€E % 



Since 

K 
k=l 

we can get by the superposition property that {Pi,i G V} is the same 
as the potential distributions {Vi,i G V} if we set the potential of the 
boundary points as V Xk = p Xk , k = 1, 2, • ■ ■ K. Because they both follow 
a harmonic function with the same boundary conditions. 

□ 



3. Proof of the main Theorem and Remarks 

In this section, we, in fact, prove a stronger result than Theorem ll.il 

Theorem 3.1. Let N = {G,c) be an electrical network with G being 
proper and C\ < c^ < C2 for all e = (i,j) G E. IfT = {x±,X2, ■ • • ,Xk} 
is boundary vertices and the boundary potential as V Xk = p Xk , < p Xk < 1 
for k = 1, 2, • • • , K, then the potential distribution of is 

V i = ^ PxkCxk +0((hm)- 1 ) 

for each i G V \ T. 

Proof. Let N' = (G', c) be the induced network obtained from iV = (G, c) 
by the induced graph G' = G[V\T] = (V',E'). Then by Lemma T2.3I and 
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equations (1), (2) and (3), we have that for 1 < k < K and i G V, 

Pi = Pr(r Xk (i) = min 1 < s < K {T Xs {i)}) 

= Pr ( T x k (i) = rnin 1 < s < K {T Xs {i)}\<ff)Pri<ff) 

+ Pr(r Xk (t) = rm ni < s < K {T X3 (t)}\^)(l - Pr{£)) 

= Pr(r Xk (i) = min x < 8 < K {T x ,(i)}\£) + Odlnn)' 1 ) (by (2)) 

= J2 Pr ( T -^ = rniniKsKKfa.tyyi^Xi^o) = j)Pr(X iiN (2t ) = j\S) 

+ 0((lnn) _1 ) (time homogeneous and Markov Property) 

= Y,Pr(T Xk (j) = TO< s <^{r X3 0')})Pr(X^(2i ) = + O^nn)- 1 ) 

= (l + 0((lnn)- 1 ))^^Pr(X^(2to)=j) + 0((lnn)- 1 ) (by (3)) 
= (l + 0((lnn)- 1 ))^^( J ) + 0((lnn)- 1 ) (by (1)) 

= Y. P ^N'{3) + 0{{\nn)- 1 ). (by P* < 1) 

if? 

Hence 

K 

Vi = Pi = Y,P* k P ? = V c + O(Onn)- 1 ). 
k=i 

Since the total current flowing into the network is equal to the current 
flowing out, we have 

K 

E E 04,-^ = 0. 

fc=l (x k ,i)eE 

Then 

if 

E E (Px.-Vc-oainnr^cto^oa), 

k=i (x k ,i)eE,igr 
which implies 

K 

E E (p* k -V c -0{{\nn)- l ))c iXk = 0{l). 

k=l (x k ,i)€E 
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Hence 

K 



k=l 



Using P5, we have 



V c = ^=lPx k Cx k + (( lnn )-l) 

This completed the proof. □ 



Let us consider a special case of theorem 13.11 in which we assign 
unit conductance for each edge and choose exactly two vertices as the 
boundary set. 

Corollary 3.2. Let N = (G, c) be an electrical network with G being 
proper and c,ij being unit conductance. If {x, y} is the boundary set and 
is added to their unit potential difference as V x = 1 and V y = 0, then the 
potential distribution of Vi is 

Vi = - f (x) „ . +0((\nn)- 1 ) 

for each i G V \ {x, y} . 



Proof of Theorem 11.11 : It is easy to see from Lemma 12.21 and Theo- 
rem 13.11 that Theorem 11.11 holds. 

Remark 1 Let us now see the concentration of potential distribution 
from a different point of view. We consider a generalized consensus model 
on G G G(n,p) with K leaders V = {xk, 1 < k < K}. For each % G V let 
Si(0) denote the score of the zth agent towards some event at the initial 
state. Set s Xk (0) — pk, 1 < k < K, at each step all agents except for 
the leaders change their scores by simply averaging the scores of their 
neighbors. Let s(t) = {s\(t), S2(t), ■ ■ ■ ,s n (t)} T be the score vector at 
step t. Then 

s(t+l) =Ps(t), 

where P is just the transition matrix of a random walk on G with T as 
absorbing states. It is known that that s(t) will approach to a vector 
s(oo) = {si(oo), 52(00), • • • , s n (oo)} T as t — >■ 00. Moreover s(oo) follows 
a harmonic function with the same boundary conditions as the potential 
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distribution on the electrical network in Theorem 13.11 when we set unit 
resistance for each edge. So the limiting score vector also has the concen- 
tration property while we choose connecting probability p large enough 
since the harmonic function has unit solution. 

Remark 2 In Theorem ll.l[ in order to make G be connected, we 
have to set the probability p larger than — . Otherwise G may be not 
connected (see [2j). But We can still consider our model on the giant 
component of G for ^ < p < ^ (see |2|). 

However while p is small, there will be many vertices on the tree tops 
which are meaningless for our model. In order to avoid this we may 
consider our model on a special case of small-world network . We can get 
a connected graph G by simply adding a random graph G(n, p) to a circle. 
In this case while p > we can still get the concentration property by 
the same method we used in theorem 3.1 since G is also proper whp. 
However while p is small we are not able to give a rigorous result now. In 
the section below we will do some simulations on small-world networks 
while the connecting probability is small. 



4. Further Discussions and Problems 

In this paper, we present the potential distributions of an electrical 
network on proper graphs and the resistance on each edge being bounds. 
It is natural to ask what the potential distributions on other graphs and 
different resistance. In this section, we consider the potential distribu- 
tions of the electrical networks on with different graphs, such as circles, 
and the small- world networks (see [19|) and the resistance Cij be i.i.d 
random variables for each (i, j) e E which t may be closed to or +oo. 
Up to now, there is no theoretical results as Theorem 11.11 since there 
seems no methods to deal with these problems. But the simulations on 
these questions may appeal some ideas. 

First, we note here if the potential distribution except for the bound- 
ary vertices are very close to a constant V c , then similarly as we proved 
in theorem ll.l[ 

V < = V c 

l^k=l ° x k 

where {xk, 1 < k < K} are boundary vertices with V Xk = p Xk . So we can 
use V c as an approximation of V c . 

We divide our simulations into three parts according to the structures 
of networks and three different independently random distributions. 



11 



Case 1: Circle. Let G be a circle on 1000 vertices, each vertex V{ has 
exactly two neighbors and V\ is connected with fiooo- Set the bound- 
ary potential as V\ = 1, V251 = 0.3, V501 = 0.7, V751 = 1. In Figure 1, 
we plot three pictures of potential distributions according to choices of 
conductance, where c^- is unit conductance, unit U(0, 1) distribution, 
and power-law distribution (seejjj), respectively. Here we use power-law 
distribution with density function as 

f(x) ~ x~ 2 - 5 for x>l. (10) 




(a) Cij is unit conductance (b) follows U(0, 1) dis- 
tribution 




(c) Cij follows power-law 
distribution 

Figure 1: Potential distribution on circles 

From Figure 1, it is easy to see that there exist no concentration 
of potential distributions on circles no matter how we choose any distri- 
bution of conductance. 

Case 2: G(n,p) model. We choose G G G(n,p) with n = 1000 , 
p = 0.01 and the expected average degree is 10. So G is proper whp. 
Set the boundary potential as V\ = 1, V251 = 0.3, V 50 i = 0.7, ^751 = 1. In 
Figure 2, we plot three pictures of potential distributions on iV = (G, c), 
where c is unit conductance, unit U(0, 1) distribution, and power-law 
distribution (see (H ) , respectively. 
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(a) Cij is unit conductance (b) Cij follows U (0, 1) dis- 
tribution 



N i Dirjt ii tMulM H ii M i i i lf ii 



(c) Cij follows power-law 
distribution 



Figure 2: Potential distribution on G(n,p) graphs 



From Figure 2, it is easy to see that concentration of potential distri- 
bution appears when we choose unit conductance as we proved in Theo- 
rem II. 1| and we can use V c as an efficient approximation of V c . Even if 
the conductance follow certain distributions such that it approaches to 
or +oo, we can still find concentration properties. But we are not able 
to give a rigorous mathematical proof in this case. 

Case 3: The small world network. We choose G to be a random 
graph G(n,p) adding to a circle of size n. Here we choose n = 1000, 
p = 0.001 so that G may not be proper. Set the boundary potential as 
V\ = l,V25i = 0.3, V 5 oi = 0.7, V751 = 1. In Figure 3, we plot three pic- 
tures of potential distributions, where c is unit conductance, unit {7(0, 1) 
distribution, and power-law distribution (see[l]), respectively. 

From Figure 3, it is easy to see that even if we choose the connecting 
probability p very small, there also exists a concentration of potential 
distribution on small-world network except for a few vertices. We guess 
this is because the random walk on small-world network also has short 
mixing time as the random walk on G(n,p) model with large p which we 
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mentioned in equation (1). 

It seems from the simulation results that only the structure of the 
electrical network will affect the concentration property. So we propose 
the following two questions: 

Problem 4.1 Does the potential distribution for an electrical network 
iV = (G, c) concentrate where G is proper and is i.i.d? 
Problem 4.2 Does the potential distribution for a electrical network 
iV = (G, c) concentrate where G is from small-world network with low 
connecting probability and cy is unit conductance or i.i.d as unit [7(0, 1) 
or power-law distributions? 



(a) Cij is unit conductance (b) cy follows U(0, 1) dis- 
tribution 




(c) Cij follows power-law 
distribution 



Figure 3: Potential distribution on small-world networks 
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